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ABSTRACT 


Solar radiation pressure exerts a mechanical force upon the surface of 
a spacecraft which intercepts the stream of photons coming from the Sun. 

For high-precision spacecraft attitude control and orbit determination, it is 
necessary to generate a precise mathematical model of the solar radiation 
force and the moment of that force; such a model must be more accurate 
than the currently used "flat surface" model, based on the radiation force 
on the effective cross-section area of the irradiated body. 

In this report the general expressions for the solar radiation force and 
torques are derived in the vectorial form for any given reflecting surface, 
provided that the reflecting characteristics of the surface, as well as the 
value of the solar constant, are known. An appropriate choice of a spacecraft- 
fixed frame of reference leads to relatively simple expressions for the solar 
radiation forces and torques in terms of the functions of the Sun- spacecraft- 
Earth angle. 

The advantage of such a model over the standard flat- surface model is 
obvious, and it is very easy to find the expressions for the error of the stan- 
dard model for any given reflecting surface. Another advantage of the model 
is that it can be used for the effects of the air drag, solar wind, etc. 
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I. INTRODUCTION 

To derive the expression for the mechanical force exerted by radiation 
upon an intercepting surface, let us consider an elementary flat surface area, 
perpendicular to the direction of the incoming radiation. Let S (Fig. 1) be 
such a surface and let J be the radiant energy per unit area and per unit of 
time, impinging on S. One part of this energy, (1 - Y)J, Y < 1, will be 
absorbed by the surface and reradiated isotropically into space. The remain- 
ing part, YJ, will be reflected according to a certain reflection law f(0), 
where 0 is the angle between the direction of the reflected particle and the 
normal to the surface area S. If the reflection is purely specular, f ( 0 ) will 
be the two-dimensional Dirac delta function. If the reflection is totally 
diffuse, the reflected particles obey Lambert's reflection law (Ref. 1): 

f(0) = cos 0. 

The total reflected radiation is in fact the combination of these two types of 
reflection. If we denote by c r the surface area of the hemisphere shown on 
Fig. 1, and by dco the elementary solid angle 

dco = sin 0 d<|> d 0, 
the total reflected radiation is 


YJ 



I dco 


where I is the radiant flux per unit solid angle on the hemisphere cr. 
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Now set 


I - l n f(0) 


where 1^ is a constant. Subsequently, integrating over the area of the 
hemisphere, we find 


VJ 


»2tt rn/Z 

••q ■ d<t>l f(0) sin 0 d0 
'0 JO 


I 0 A(f), 


where 

( 1 for specular reflection 
w for diffuse reflection. 

The momentum of the radiation is J/c, where c is the speed of light. 
The momentum exchange due to the reflected beam of photons in the direction 
of the normal to the surface area S is 




f(0) sin 0 cos 0 d<f> d0, 


and the tangential component of the momentum is zero. After an integration 
over the surface of the hemisphere <r, the integral above becomes 



— cos 0 dw = 
c 



sin 0 cos 0 d0 



where 

! 1 for specular reflection 
2/3 for diffuse reflection. 
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The total pressure is the sum of momentum exchanges due to the incident 
and reflected radiation 


P = 



Since J is inversely proportional to rgp (fgp being the heliocentric 
distance of the irradiated body) the acceleration due to the radiation pres- 


sure is 


a SRP = 



where 


k 


£p <, 1 + VB(f) 
c M 


Here, Jq is the value of J at the distance of the Earth from the Sun (one 
astronomical unit), known as the solar constant: 


3 2 

Jq = 1. 353 X 10 watts/meter (Ref. 2), 


while M is the mass of the reflecting body and S is the irradiated surface 
area. 
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Fig. 1. Radiation reflected from an elementary surface area 
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II. SOLAR RADIATION FG*vCE ON AN INCLINED SURFACE 

Let us assume now that the direction of the incoming radiation is 
inclined by an angle 0 to the normal to the elementary surface area S (Fig. 2). 
The incident radiation produces the force 

F I = f n cos 6 


where F^ is the radiation force for 0 = 0, i. e. , in the case when the incident 
radiation is perpendicular to the intercepting surface, considered in the pre- 
ceding section. This force is given by 


F 


N 


k srp_ 

__ s 

r SP 


( 1 ) 


where K OOT1 is the solar radiation pressure constant. The value of this con- 
SRr 

stant can be easily calculated from the solar constant Jq- It is 

„ < AU)2j 0 „ . 2 

k srp ' — 5 kg m/s • 

Taking for the value of the astronomical unit (Ref. 3) 

AU = 149,597,893.0 km 
and for the speed of light (Ref. 3) 

c= 299,792.5 km/s 


we finally compute 


K SRp = 1.0227 X 10 17 kg o/s 2 . 


( 2 ) 
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The component of the force Fj along the normal to the surface area S is the 
pressure force (Fig. 2) 


Pj = Fj cos 0 = F^ cos 2 0, (3) 

and the component of the same force in the tangential plane to the surface is 

Tj = Fj sin 0 = Fj^ sin 0 cos 0. (4) 

As mentioned beforehand, of all the photons impinging on the surface S, 
the portion of 1 - v remains absorbed and reradiated isotropically into the 
neighboring space, while the remaining portion Y is reflected specularly or 
diffusely. Let py be the portion of the specularly reflected photons. The 
force produced by the reflection, after the collision with the surface S, is 

f r = ewy 

and the components of this force along the normal to the surface of collision 
and in the tangential plane are 


P R = f r co, e | 

T R = F R sin 8 ) 

so that, substituting F R and Fj, we finally find, in terms of F^, 

P R = PYF n cos 2 0 (5) 


and 


T r = PYFjy sin 0 cos 0. 


( 6 ) 
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The force exerted on the surface S by the diffusely reflected photons, (1 - (3)y, 
is perpendicular to S. Its magnitude is 


2Y 

F D = -<i-P>f n «. e. 


(7) 


The total radiation force along the normal to the surface is then the sum 


N 


= P I + 


P R + F 


D’ 


and the total tangential force is 


T = 



Substituting P T , P_, F_^ T t , and T from Eqs. (3-7), we obtain, in terms 
1 r\ L) 1 

o f F N . 

P N = F N [(1 + pY) cos 2 e + -f^ 1 - P) cos e ] (8) 

and 

T = F n (1 - pY) sin 6 cos 6. (9) 

The polygon of forces and the resultant force R are shown on Fig. 2. The 
magnitude of the resultant force is 

r n * f n cos 6 

where 

2 2 4Y^ 

Kj = 1 + pV +^-(1 - 

K 2 = - Y){1 + pY) 

K 3 = 2PY- 


^Kj + K 2 cos 6 + K 3 cos 26 


( 10 ) 


P) 
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The ang;e between the direction of the resultant force and the normal to the 
surface (oriented toward the surface) is 


6 


R 


arctan 



(1 - PV) sin 6 
- P) + (1 + pY) cos 6 


(ID 


8 
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Fig. 2. Polygon of forces acting on an elementary surface area 
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HI. COMPUTATION OF THE SOLAR RADIATION FORCE 


Consider an infinitesimally small element dS of the irradiated 
reflecting surface S (Fig. 3), and let N* be the unit vector of the normal to the 
surface at the point A (x, y, z) oriented in such a manner that 0 < 6 S 90 deg. 

If 4>(x.y,z) = 0 is the reflecting surface, then, obviously, 


N* 


grad $ 
|grad ®|* 


and the oriented vector element of surface is 


( 12 ) 


dS = N* dS = e^ dy dz + dz dx + e^ dx dy, 


(13) 


where ej, e^, e^ are the unit vectors of a reference frame xyz. Hence it is 
easy to see that 


dy dz _ dz dx _ dx dy 

TT* • « 1 TT* • ej T7* • e 3 


(14) 


The unit vector of the direction of the tangential force T* lies in the 
plane containing the normal, incident and reflected rays, which pass through 
the point A. Also T* is directed toward the half of the plane in which the 
source of the light is, so that the angle between and the direction of the 
light source is acute. Denoting by u the unit vector of the Spacecraft-Sun 
direction, we find, from the conditions 

T* • N# = 0 
(uXN*) • T* = 0 

u • T* = sin 0 


that T 7 * is given by 


T* = -TTf - N* cot 0. 
sin 6 


(15) 


10 


JPL Technical Memorandum 33-494 



Denote by dF the resultant radiation force acting upon the element 
of surface dS. Since 


dF = -N* dP N - T* dT, 


and since, from Eqs. (8) and (9), 


dP N = [(1 + pY) cos 2 6 +-y(l - p) cos eJdFj 


dT = (1 - pY) sin 6 cos 0 dF. 


where, from Eq. (1), 


SRP 


K(r gp ) dS, 


we can write, for the elementary radiation force, 


dF = -K(r sp )|p(e)N* + t(0)T*]dS. 


The total solar radiation force over the whole area S is then 


F = -K(r sp ) | / [p(0)N* + t(0)T* ] dS, 


p(0) = (1 + PY) cos 2 0 + ~ (1 - P) cos 0 
t(0) = (1 - pY) sin 0 cos 0. 
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The element of surface dS is given by Eq. (14), and the choice of the 
projection of the area S on the coordinate planes, over which the double 
integration should be performed, is completely arbitrary. One should 
choose, of course, the coordinate plane in which the projection of the area 
S is the simplest geometrical figure which renders the technically simplest 
integration procedure. 

Replacing, in Eq. (18), vector T* by its value in Eq. (15), we give 
the expression for the total radiation force the following form 



F = -K //|B(e)N*+^|L- 


■] 


dS 


where 


B(e) = p(e) - t(e) cot e, k = K(r gp ). 


Using the notation 


2Y 

Cj = xd - 13) 
C 2 = 1 - PY 
C 3 = pY 


we can write the function B(0) in the form 

2 

B(0) = Cj cos 0 + 2C^ cos 0, 


( 20 ) 


( 21 ) 


and the radiation force becomes 



( 22 ) 


12 


JPL Technical Memorandum 33-494 



The second integral in Eq. (22), because of 


cos 0 = u • N*, 


becomes 



//■ 


(u * "FT*) dS 



u • S. 


Denote by A , A , A the projections of the total illuminated area 
7 yz xz xy r J 

on the coordinate planes. Subsequently, 


S = 


A e. 
yz 1 


A e~ 
xz 2 


+ A e- . 
xy 3 


If a is the angle between the direction of the incident Sun ray, defined by thr 
unit vector u, and the z-axis of the system, then 

u = e^ sin a + e^ cos a (2°) 


and 


u*S = A sin a + A cos a, 
xz xy 


The total radiation force is then 


jj. 


F = -K I |B(e) dS -KC,(A sin a + A cos a) u 
' 2' xz xy 


(24) 


(25) 
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where, from Eq. (21) 


B(e) = •— (u • N*) |l - p + 3 P(xi • N*)J, 

and also where 


(26) 


N* = 


grad $ _ _ V t 

jgrad "*j ~ |V$| " W ’ 


(27) 


since $(x, y, z) = 0 is the equation of the illuminated surface area and 
W = IV*|. 


Incidentally, for a black body, y = 0, (3 = 0, C ^ - 1. and the radiation 
force (Eq. 25), is given by 


F = -K(A sin a + A cos a) u. 
xz xy 


(28) 


The first integral in Eq. (25) in the form 



may be hard to evaluate since all three projections of the area S do not 
necessarily have to be simple geometric figures. Therefore, as wa3 
mentioned before, the best method for the practical evaluation of the double 
integral in Eq. (25) would be to use it in the form 




B(e)N* dS 


s 
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and to choose di> from one of the three expressions given by Eq. (14). Hence 
we can adopt for the total form of the total radiation pressure, the following 
expression: 


F 



1 - P + 3|3 


(u • V $) 
W 


( u • vs> ) 

w 2 


V$ dS 


- (1 - PYJK^-— ^(A xz sin a + A xy cos a)u 
where, from Eq . (14), 


or 


dS 


W dy dz 
V$ • | 


W dx dz 

R* • e 2 | 


W dx dy 
I V4> - e 3 J 



— dx dy. 

II# 

| a ^ 


(29) 


(30) 


In the practical application of Eq. (29), one should be aware of the fact 
that the normal to the surface, defined by the unit vector N* f bisects the 
angle between the incident and the reflected solar rays so that 0 < 0 < 
and, therefore, 


N* • u > 0 


or 


VS> *u > 0. 
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Hence, the direction of should be chosen in such a manner that the 
inequality above is aiways satisfied. At the end of this section, it should be 
noted that from cos 0 = u • N* and Eq. (27), 


cos 6 
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Fig. 3. Orientation of unit vectors along 
the tangent and normal to the 
reflecting surface 
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IV. MOMENT OF THE SOLAR RADIATION FORCE 

The same considerations used in the preceding section can be extended 
to the problem of finding the moment of the solar radiation force via a double 
integration over the illuminated surface S. The elementary moment of the 
force dF, acting upon the elementary surface area dS, is given by 

dM = 7 X d F 

where r is the position vector of the point A (Fig. 3) in the system of 
reference xyz. From the expression for the elementary force in the form 


dF = -K ^B(0)N* + C 2 u cos ej 


v/e find that 


dM = -K [B(0)(r X N*) + C 2 (r x u) cos ej dS, 


so that the total moment of the solar radiation force over the area S is 


//■ 


M = -K I J B(0)(r X dS) + KC 2 u X I I r cos 6 dS 


ft 


M = -K 


f/ B(9,< r X dS) + KC 2 u dS ) 
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Finally, setting 


ds = at dS> 


we obtain 




1 - p + 3p 


(u 


7, dS 

w J w z 


+ (1 - 0Y)K 




dS. 


(31) 


The unit vector u is directed toward the Sun. Denoting the heliocentric 
position vector of the spacecraft (Sun-Probe vector) by r CD as before, we 

Or* 

have 


u 



(32) 


Equation (31) is the final form of the expression for the moment of the solar 
radiation force. The quantity A is, as before, the projection of the illumi- 
nated surface area on the xy coordinate plane. 

Before moving to the next section, in which we shall proceed with the 
integration of Eqs. (29) and (31) for the spacecraft components, we must 
mention that the reflection does not strictly obey Lambert's law for some 
materials and, therefore, the function f(0) should be determined in such a 
manner that it represents the actual response of a certain material to 
reflection. The parameter B(f), mentioned beforehand, can then be com- 
puted using f(0) for that particular material. 

Another factor, which can also contribute an extra force of small 
magnitude, is the reradiation of the thermal energy. In Ref. 4 the author 
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suggests an alternative expression for Eq. (8), which accounts for the 
above-mentioned effects* The expression is 

P N = F N j(l + PY) cos 2 e + B(f)[(l - p)Y + k(l - Y)J cose 

where is given by Eq. (1) and where the parameter k depends on the 
temperatures and emissivities of the front and back sides of the reflecting 
surface and on the absorbed and conducted flux. For adiabatic surfaces, 
k = i. 

It is obvious that with the expression above for P^, the function B(6) 
in Eqs. (29) and (31) will have a somewhat different form. Namely, because 
in this case p(9) is given by 

p(e) = (1 + J3Y) cos 2 e + B(f) [(1 - P)Y + k(l - Y)J cos 6, 
while t(0) remains unchanged, and 

b(g) = p(e) - t(e) cot e. 


w r e find that 


Bjte) = 2/3Y cos 2 e + B(f) £(1 - P)Y + k(l - Y)J cos 6. 
Hence, we see that, using the notation 

c\ = B(f) [(1 - J3)Y + k(l - Y)J, 

we can write the function B(0) in the form 


BjO) = c\ cos 0 + 2pY cos 2 0. 


20 
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Comparison of this expression with B(6) given by Eq. (21) indicates 
that B j(0) and B(6) are one and the same function of 0, differing only in the 
value of the constant Cj. This difference is 

c\ - Cj = [b( 0 - |] Y(1 - p) + kB(f)(l - Y). 
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V. THE SOLAR RADIATION FORCE ON THE REFLECTING SURFACE 
OF A PARABOLIC ANTENNA REFLECTOR 

The high-gain antenna of the Pioneer spacecraft has the shape of a 
paraboloid of revolution, with its axis of symmetry presumably pointing toward 
the Earth. The configuration of the spacecraft and the position of the antenna 
reflector relative to the other components is shown on Fig. 4. A simplified 
sketch of the reflector only is given on Fig. 5. , displaying the dimensions of 
the dish: the semidiameter 6 and the depth C. 

Since the double integration, indicated in Eq. (29), constitutes the major 
technical difficulty in the derivation of the expression for the solar radiation 
force, it is a matter of utmost importance that the coordinate system, with 
respect to which the integration should be performed, be properly and ade- 
quately chosen. 

The most suitable frame of reference, which would be most objective 
for the purpose of the double integration, provided that the roll-axis of sym- 
metry of the spacecraft coincides with the spacecraft - Earth direction, can 
be defined in the following manner. Taking the z'-axis along the spacecraft — 
Earth direction, we find that the unit vector of the z'-axis is given by 



where Fgp is the geocentric position vector of the spacecraft and = 

I f epI- 

The y'-axis of this system lies in the plane defined by the unit vector 
k and the unit vector u, defined by Eq. (32). The x'-axis is perpendicular to 
this plane; hence 

t u X E 

1 = |iTxT| 

and 

j = k X 7. 


22 
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Denote by the angle between the positive direction of the z'-axis and the 
direction of the unit vector u (Earth- Probe -Sun angle). Then 

|u X E| = sin 


Taking the expression for u from Eq. (32), we can write explicitly 


1 = 


r SP X r EP r SP X r EP 


r SP r EP Sin Q 0 | r SP X r EPl 


7 r EP X (r EP X r SP ) 
J " 2 

r EP r SP Sin “o 


k =- 


EP 

r EP 


This frame of reference is noninertial since the two vectors Tgp and 
vary with time. 

With the introduction of the heliocentric position vector of the Earth, 


R E r SP " r EP ’ 


whence 


r EP " r SP ‘ R E ' 


the formulas above can be simplified and reduced to the form 

T- r e x f 

1 rR E sin ^ 
j = k X i 

— — _ T 2 2 -i-l/2 

k = (R e - r)Jr + R e - 2rR £ cos 4 J 
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where r = r gp» r = |^|» and 4* ** the angle between the vectors r and Rj, 
(Earth-Sun- Probe angle). 

Denote by aj, a^ the unit vectors of direction of coordinate axes of 
an inertial reference frame (for example, the 1950 equatorial coordinate sys- 
tem), and let X, Y, Z, and Xg, Yg, Zj, be the components of vectors 7 and 
Rg, respectively. The instantaneous angular velocity vector w, defining the 
rotation of the above-described rotating frame relative to the inertial equa- 
torial frame of reference, is then 


u 




+ k X 



The unit vectors of the noninertial system of reference can easily be 
expressed in terms of the coordinates of the spacecraft and the Earth and 
their derivatives in the inertial frame. Thus 




24 
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where 


p = [( X - x e) 2 + ( Y - Y e) 2 + < Z - Z E> J 


1/2 


r = (X 2 + Y 2 + Z 2 ) 


1/2 


-(■ 


R E = (4 + 4 + Z E 


, 1/2 


cos 


4 » = 


xx E + yy £ + zz E 
7SZ 


However, due to the existence of torques, the axis of symmetry of the 
reflector moves relative to the above-described system of coordinates, 
defined by unit vectors i, j, k. Therefore, for the computation of the radia- 
tion force from Eq. (29), by double integration, we can establish another 
coordinate system, defined in the following manner: 

(1) The z-axis is taken along the axis of symmetry of the reflector. 
The unit vector of this axis is *e^. 

(2) The y-axis, defined by the unit vector e^, lies in the plane of 
vectors e^ and u. 

(3) The x-axis is perpendicular to the yz-plane, forming a right- 
hand-oriented triad with the first two axes. Its unit vector will 
be denoted by"ej. Therefore, 

_ 1 X-g 3 

e l sin o 


where a = ^.(u, e"_ ), namely, the angle between the z-axis and the direction 
of the incoming radiation. The xyz-coordinate system described above is 
shown on Figs. 5 and 6. 
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The geometry of the parabolic reflector is shown on Fig. 6. Placii g 
the origin of the above-described coordinate system at the vertex of the 
paraboloid, we can write for the equation of the reflecting surface 

4>(x, y, z) = z - X.(x 2 + y u ) n 0 (33) 

where X is the parameter 

^ ~ , (34) 

6 ^ 


depending on the dimensions of the reflector only. For the Pioneer F/G 
spacecraft, we have, approximately, 

X = 0,202 meters"^ . 

The magnitude of the illuminated inside area of the reflector is a func- 
tion of the angle a . For small values of a, the whole area is illuminated and 
the projection of this area on the xy-plane is a circle of radius 6 . When the 
solar ray passing through the point A on the brim of the reflector dish (Fig. 
6) becomes tangent to the surface, the tangent line cuts off the segment C on 
the negative z-axis, so that 


a = 90° - ft 

whe re 

tan fi = ~ 

6 


or 


where 


tan ft = 2 tan , 


tan \\i = 


£ 

6 * 


The two angles, ft and are shown on Fig. 5. 


(35) 


(36) 


(37) 
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If a > 90° - ft, one part of the inside area of the reflector will be 
in the shade, and the projection area in the xy-plane will not be the full area 
of a circle, but a somewhat more complicated figure. To find this figure, 
let us start with the equation of a solar ray passing through any point (x^, 
y Q , Zq = £) on the brim of the dish: 


where 


x ^ x o = i : y o = z - 4 

0 sin a cos a 


2 

x 0 




Take y^ = 6w, where w is a parameter. Then, in terms of u>, the equations 
above yield 



y = w6 + t sin a 
z = C + t cos a . 


(38) 


These equations represent, in the parametric form, the equation of the 
inclined elliptic cylinder, shown on Fig. 6. The semimajor axis of this 
cylinck _ is 6, and its semiminor axis is 6 cos a. The line of intersection of 
this cylinder with the paraboloid of revolution 


z = \(x^ + y 2 ) 


occurs for the values of t which are the solutions of the quadratic equation 

(X. sin^a)t^ + (2X6w sin a - cos a) t + (\o^ - £) = 0 . (39) 
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2 r 2 

However, from Eq. (34), - b =0, and the solutions of Eq* (39) 


are 


h * *2 


=. 26 w) . 

smo'\ a, / 


For tj = 0, we have 


= 6\lT 


w 


y = 6w 


and, eliminating the parameter w, we obtain the equation of the circle 


2 2 .2 

x + y = o . 


For 


t. = 


sin a 


(cot a cot 4* - 2w) , 


we obtain the equation of the line of intersection in the parametric form. 


x = 6 V 1 - w 2 


y = 6 (cot a cot - w) 


(40) 


= 6 tan \!f> J(1 - w^ - ) + (cot a cot 4» - w)^J 


The projection of the curve given by Eq. (40) on the xy- plane is given 
in the parametric form by 


x = 6 \ 1 - 

y = 6 (cot a cot ^ - w)# 
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Eliminating w, we find 


(41) 

(42) 

The two circles intercept if 0 < y^. < 26 or, in other words, if 

0 < cot a cot »|* < 2 , 


x 2 + 


(y 


.2 * 2 

y c } = ^ 


where 




= 6 cot a cot ^ 


or 


or if 


0 < cot or < 2 tan 4* = tan £2 . 


90° - £2 < a < 90*. 


For a < 90° - £2, the whole inside area of the reflector is illuminated. 
For a = 90*, the whole area is in the shade. The y-coord ; nate of the 
point B (Fig. 6) is 


y B = y n - 6 = 6(cot a cot 4> - 1) , 


(43) 


and the two circles intercept in the points with coordinates 


= ±2 -^4 - cot 2 a cot 2 4* 


yg = 2 cot a cot 4* 


(44) 


The shaded area on Fig. 7 represents the projection of the illuminated 
inside area of the reflector on the xy- plane. The angle is given by 


y s j Aq 

! n 4 > q = -jr- = — cot a cot 4* = = cot a cot £2 
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where 


Ag = cot a cot 4 1* = 2 sin 6 q 


2 2 2 

With this denotation, the equations of the two circles, given by x + y =6 
and Eq. (41 ), are 


y = ±yji> 2 - x 2 


y = 6 Aq ± yj S 2 - x 2 


and the parametric equation of the line of intersection (Eqs. 40) becomes 


x = 6 M 1 - w 


y = 6(A n - w) 


= ^[(1 - w 2 ) + (A q - w) 2 ] I 

= i(l + Aq - 2A Q w) . J 


The projection of this curve on the xz-plane is, in the parametric form. 


x = Wl - w 


z = til + Aj - 2A 


o w ) • 


and, after the elimination of the parameter w, we obtain the ellipse 


Z = + Aq - - X 2 
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or 


= 4(1 + Aq) - 2 cot atyj 6^ - , 


or 


(f ♦ 


, 


26 cot « 


( 48 ) 


The projection of the same line on the yz-plane is, in parametric form, 


y = 6 (a q - w) 


= + A o ■ 2A o w ) 


and, after the elimination of the parameter w, we obtain the equation of the 
straight line 


= 2y cot a + 4 |l - j . 


(49) 


For z = 4, we get 


6A q 

y = — = y s 


and for the point of intersection with the parabola z = Xy (point B), from 


we find 


Xy^ - 2y cot a - 4^1 - A^j = 


y = 6(A 0 - 1) = y B 


which is the value already obtained in Eq. (43). 
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Fig. 


The projections defined by Eqs. (48) and (49) are shown on Fig. 8 and 
9, respectively. 

From the equation of the reflecting surface, we find 


which yields 


and 


so that 


2 2 

S(x, y, z) = z - \(x + y ) = 


VS = -2\xej - Z\ye^ + e^ 


W 


= I VS | = 


1 + 


4\ 2 (x 2 


2 . 

y ) 


N* = — (-2X.xej - 2Xye 2 + 


e 3> 


From 


cos 0 = N* • u 


where u is given by Eq. (23), we find 


_ cos a - Z\y sin a 
cos 0 = yj- 1 


(50) 


(51) 


(52) 


(53) 


where 0 is the angle between the incident solar ray and the normal to the 
surface $ at a particular point. Therefore 

cos 0 > 0 


or 


cos a - 2Xy sin a > 0 . 


This condition is satisfied for 


. cot q 

y < — 



* y s- 


which is correct. 
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Returning to Eq. (22), we shall first find the force over the whole area 
of the reflector,- when a < 90® - S2, that is, when the entire interior of the 
mirror is illuminated. Introducing the polar coordinates 

r r r 

x ~ 2\ cos 4*’ Y = 2\ sin dx dy = — g dr d«(> , (54) 

4X. 

we can write 


0 _ cos a - r sin <)> sin at 

COS U — "• 


yjl + r' 


(55) 


and 


Since 


dx dv r dr <l<t> L . T 

ds = iijrT Fjl = “^2 V 1 + r 


dS = N* dS = f dr , 

4\ C 


(56) 


Eq. (22) becomes 


F = - 


K 

4X 2 


ff KC 2 — ff I 1 

JJ™ B(0)r dr d<j. - —j u JJy] 1 + r cos 0 r 


dr d<(> . 


xy 


xy 


Hence, 


F = 
x 


K 

— 2 
4\ 


K 


// 


B(0) r^ cos 4> dr d<|> 


xy 


F = — A 

y 4 k 


fj me) t z si 


(57a) 


in <)> dr d«|> - KC 2 I sin a 


xy 
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F = - 
z 


B(0)r dr d<J> - KC^I cos a 


where 


(57b) 


1 + r cos 0 r dr d(f> 


with B(0) given by Eq. (21). 

The integration of the right-hand sides of the three equations given 
above is rather complicated in polar coordinates, particularly when y^. < 6\/2, 
or when (Fig. 7) 


<J> 0 < 45°, 


because the radius vector OS^ (or C^) intersects the upper circle on Fig. 7 
in two points. The integration is simpler in rectangular coordinates, as will 
be shown later. 

However, before going deeper into the details of the integration over 

the area A , we shall first derive the components of the radiation force 

over the whole inside (front) area of the parabolic reflector, when a < 90° 

- ft. The projection of this area on the xy- plane is the area of the circle 
2 2 2 0 

x + y = 5 , which we shall call A . The components of the force of radia- 

tion over the whole area will be denoted by 


F° = {f°,F°.F 0 }. 

» x* y’ z J 


The limits for the integration over the area A ^ are: for 4>» <f>j = 
= 2ir , and for r, r^ = 0, r~ = 2X.6 = 2 tan ^ = 2 m. Then 


34 


JPL Technical Memorandum 33-494 



<■5 IT 


(cos a - r sin <j> sin a?) 2 — — cos <J> dcj> 

1 + r 


4\ 


» 17 


KC, I -2 ’ r 2m p 2 dp 

(cos a - r sin 4> sin a) cos <j> d<J> 




1 + r 


( 58 ) 


2ir /’2m 


F° = 

y 


-3.fi 


(cos a - r sin <f> sin a) 2 — — sin <j> d<j> 

1 + r 


4\ 


» a: 


KC, /‘ 2ir f 2 ™ r 2 dr 

(cos a - r sin <f> sin a) = - =^ sin <|> d<^> - 'KC^l sin a 


'JT7? 


(59) 


0 kc 3 r 2 " r 2m 2 p dp 

* 77 Jo Jo 777 


KC 


4\ 


n 2m 


r dr 


(cos a - r sin 4* sin a)—j==s d<j> _ KC_I cos a- 

VI * ■ * 


+ r 


It is easy to see that, from Eq. (58), 


(60) 



(61) 
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Equations (59) and (60) can be written in the form 


y = k ( C 3 I 2 i + c i I 22 - C 2 Isin «) 


= -K(C 3 I 31 + 0^22 + C 2 I cos a) 


( 62 ) 


where, after the integration is performed, 


2tt 

T sin 2a I . 2 , , , 

hi - — rr I sm * d * 


22 


“I 

^ f 
^ Jo 


sin d4> 


Jo I- 1 " 
Jo VTT7 = 


2 2 
= sin 2<*(1 + 2 cot ft In cos ft), 


(63) 


tt 6 . , 2 ft /0 , 

-—j— sin a tan —(2 + sec ft). 


(64) 


2 r Zv r 2m , .2 f 

"--zrj o Vo ~ 2 ~J 0 


*2n /»2m , 

. 2 . , . / r dr 

sin <}> d<)> 

1 + r‘ 


ir6“jjsin 2 a + 2 cot 2 £2(1 - 3 cos 2 a) In cos fij 


■ *£[« - 


2 cot 2 £2 In cos £2) - (1 + 6 cot 2 £2 In cos £2) cos 2arJ> (65) 


/ *2tt /*2m 

cos a I , . / r dr _ ,2 cos £2 

~~J I d * 777^ = 2Tr6 cos °l + coslT (66) 

4V Jo Jo ^ 1 + r 


32 ' 4\ 2 


J /*2ir /*2m 

d<0 

0 Jo 


I = 


cos 

4\ 


r dr = tt6 cos a . 


(67) 
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Collecting all terms given by Eqs. (63-67) together, we can write Eqs. 
(62) in the following form: 




z 


■ tt6 2 Ki 


|fj sin a + f^ sin Zaj 


-“rr6 2 K|gQ + g^ cos a + cos Zaj 


where 


, 2Y .. „. l - cos £2, „ . o\ 

f l - 9 ( J - P)j + cos ^( 2 + sec n )' 


f_ = 1 + 2|3Y cot 2 £2 In cos £2, 


1 2 

go = 2 “ cot ^ In cos ^ * 


4Y,, cos £2 

§1 = T (1 - P } 1 rcos’Q 


^ - pV (1 + 3 cot 2 £2 In cos £2) 


( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 


For the Pioneer spacecraft 
6 = 1. 3716m 
£ = 0.3803 m 

it 6 2 = 5. 9102 m 2 
m = tan = 0. 277267 
ip = 15*. 498 


JPL Technical Memorandum 33-494 


37 



tan ft = 0. 554535 


ft = 2 9 0 1 2 

and 

fj = 0.046754V(1 - p) 
f 2 = 0.5 - 0. 371 922fJY 
g Q = 0. 5 + 0.43596ipV 
g, = 0. 622046v(l - P) 
g 2 = 0. 5 + 0. 307883(3y ? 

To compute the solar radiation force exerted upon the inside surface 
of the reflector using Eq. (22) or Eq. (29), when only one part of the surface 
is illuminated, one must select the proper projection of the surface on one 
of the coordinate planes which yields the simplest double integration. 

From the preceding considerations, it follows that the projection in 

the yz-plane (Fig. 9) is the simplest geometrical figure composed of a part 

2 

of the parabola z = \y , A' 0B, and the segment BS of the straight line 

z = 2y cot a + 4(1 - A^) . 

Since, from Eq. (14), 

dS = n wfr | = dv dz = ^ ^ — \[\ dv dz . 

|fl * 211 


Eq. (22) becomes 
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where 2A^ the projection of the surface rea S on the yz-plane (the area 
projects twice). The equation above is equivalent to 


F 




B(6) 






dy dz 


KC- 



W cos 0 


dy dz 



The component of this force along the x-axis, because of the axial 
symmetry of the reflector, vanishes. The other two components are 


F 

y 


F 

z 


where 


2Kn/I 




. cos 0 + 2C., cos^ fi) y d y — =Jj sin a 

1 3 / I ~Z ^ 

y z - Xy 


cos 0 + 2C-, cos^ 0^ — J cos a 
3 / f ^ \/X 


\ z - \y 


(74) 


J = 



(cos a - 2\y sin a) 


dy dz 
>/z - \y 2 


(75) 


Integration of Eqs. (74) and (75) cannot be performed without diffi- 
culties taking z = z(y) because for the limit z = ^y^ the denominators of 
the integrands in Eqs. (74) and (75) vanish. Therefore, the integration 
should be performed taking y as a function of z. Thus, the limits for inte- 
gration over the part A’B^BS are 


Yl (z) = Vf 

y 2 w ■ ^ - (> - A o) c ] 
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z 


2 


= t 


and, for the integration over the remaining part B^OB, the limits are 
y^(z) = - 

y'z {z) = 


z 


1 


= 0 


z 


2 


z 


B ‘ 


Taking cos 0 from Eq. (53), we find 


KC, 


F y = + 2C 3 J 22 ) - -rfl sin « 


(76) 


K KC 2 

F z = -^ (C 1 J 31 + 2C 3 J 32> - “7T J COS « 


where 


21 


// 


yz 


CO se y - dy - d l , 
^ z - ^y 2 


22 


// 


cos 8 


yz 


2 fi y dy dz 
>Jz - \y 2 
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The evaluation of these four integrals, together with the fifth integral 
J with cos 0 taken from Eq. (53), is somewhat complicated; it is relatively 
simple to reduce the double integration to the single integration with constant 
limits Zj and z^, or z[ and zJ,* However, it may be much simpler and more 
suitable to apply directly the numerical double integration procedure with 
variable limits. 

The integration for 90 ° - ft < a < 90 ° can be performed in rectangular 
coordinates in the following way. Again taking the projection of the illumi- 
nated surface of the ieflector on the xy-plane, we can write (see Fig. 7) 




Y 2 <*> 

j(x) 


and, because of the 


symmetry of the reflector with respect to the yz- plane. 


and 



r 

L 


> n 2 w 

A,w 
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so that 



In order to avoid the double integration, an approximative method can 
be introduced. The method is based on the assumption that the ratio of the 
components of the force over the illuminated area and the components of the 
force over the whole inside area of the reflector are the same as the ratios 
of the respective projection areas. Namely, we can assume that 


f = n(<*)F° 

y y 

F = n(«)F° 
z z 


where, from Fig. 7, 


n(«) 


P rojection of the illuminated surface 
Area of the circle 


(77) 


Explicitly, this ratio is 

2$ + sin 2<J>. 

nM - 0 , 

where <|>q is a function of the angle a, defined by 

sin <|>g = i cot cot q = cot cot a . 

For the Pioneer spacecraft, = 15*498 and 

sin 4> 0 = 1 • 803 cot a . 


(78) 


(79) 
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For a = 90® - ft, when the whole inside area of the reflector is illuminated, 
sin = 1, = w /2 and q(a) = 1. Therefore, 


n(«) 


1 ^(2$q + sin 2 $q), for 90® - O < a < 90® 
1, for a < 90® - 


The graph of this function is shown on Fig. 10. 
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ASTEROID/METEOROID 
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Fig* 4. Configuration of components of the Pioneer spacecraft 





t 



Fig. 5. High-gain antenna reflector of the Pioneer spacecraft 
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Fig. 7. Projection of the illuminated inside area of the 
reflector on the xy-plane of reference 
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w 0 # 10* 20° 30° 40“ 50“ 60“ 70“ 80“ 90“ 

a 

Fig. 10. Graph of the blocking function for the 
concave side of the parabolic 
reflector 
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VI. THE MOMENT OF THE SOLAR RADIATION FORCE ON THE 
REFLECTING SURFACE OF THE PARABOLIC REFLECTOR 

Equation (31) gives the final form of the expression for the moment of 
the solar radiation force over a surface area S. For the inside area of the 
parabolic reflector, 

r K dS = r X N* = (y®^ ' xe^Hl + 2\z) dx dy 

u X r = ej(z sin a - y cos nt) + e^ x cos a - x sin a 


and the three components of the moment vector M are 


M = 
x 


-K j j'yBWU + 2\z) dx dy + KC JI" z sin a - y cos or) cos 8 dx dy\ 


xy 


xy 


■//■ 


M = K| IxB(G)(l + 2\z) dx dy + KC cos a I IxW cos 0 dx dy 


“// 


(80) 


xy 


xy 


“If 


M z = -KC^ sin a J J xW cos 0 dx dy 


xy 


/ 


To integrate over the whole inside area of the reflector, we must sub- 
stitute B(0) from Eq. (21) and cos 8 from Eq. (55) and use the polar coordinates 
defined by Eq. (54). Also, because of z = r^/4\, 


1 + 2\z 


2 + r 


2 


2 
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and 


W = v/l + r 2 • 


Since cos 0 contains sin 4> only, and B (0 ) is a function of cos 0 only, all 
integrals containing x = (r/2X.) cos 4> (the two integrals in and one in M^) 
vanish when integrated from <}> = 0 to <f> = 2ir, because 


n2m 

1 

= 0 


R(r)F(sin <t>) cos <|> dr d4> 


0. 


therefore. 




0. 


( 81 ) 


The M x component has the form 


,.0 K 

M = 7 

x 16 K 3 


c rz* rz m 3 

—y sin a I dt(> I — - — ^ 

Jo Jo \Tl + r 2 

sin 2a I 

Jo 


dr 


2 " r 2m 3 (? . 2 . 

d4.| r - (2 1 0 dr 


Jo 1 + 


•2t r /*2m 

+ sin 2a/ d$ | r dr 


JPL Technical Memorandum 33-494 


51 



or, after the integration, 


0 ttKCj 5 

M = r sino (4 + sec ft - 5 sec ft) 

x 80X . 5 


^KC 4 

+ *• sin 2 a (sec ft - 1 + 4 In cos ft) 

32 ^ 


7tKC 2 4 

+ r sin 2a tan ft. 

32\ 5 


Finally, substituting \ = (tan »4> ) / 6 , we find 


= -n-6^^K(hj sin a + sin 2a) 


(82) 


where 


hj = y^(l - P) cot^ft(4 + sec^ft - 5 sec ft)' 


h 2 = + 2f3Y cot^ft (tan^ft + 2 In cos ft)J 


(83) 


For the Pioneer spacecraft, these values are 


h 2 = 0. 334912 Y( I - p) 


h 2 = 0. 5 + 0. 4165020 


and 


tt 6 2 ; = 2. 2477 m 3 . 


52 


JPL Technical Memorandum 33-494 



For 90° - O < a < 90°, when one part of the inside area of the reflector 
is in the shadow, 



M = 

y 

M = 


q(<*) 


0 


0 


( 84 ) 


where the "blocking function" r((a) was given in Eq. (78). 
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VII. THE SOLAR RADIATION FORCE ON THE BACK SIDE 
OF THE PARABOLIC REFLECTOR 

Assuming that the thickness of the reflector is negligible, we can write 
the equation of the outside surface in the form 

$ s z - X(x^ + y^) = 0 . 

The normal to this surface is now pointing outwards, i. e. , the grad $ is now 
given by 


grad $ = 2Xxej + ZXye^ - 
and the unit normal vector N* is 


N* = ^(2Xxe^ + 2Xye^ - e^) 


(85) 


( 86 ) 


where W is given by Eq. (51). 

The geometry of the outside (back side) area of the reflector is shown 
on Fig. 11. The curve which separates the illuminated part of the surface 
from the part in the shadow is the locus of points at which the solar rays are 
tangential to the surface of the paraboloid. If 

x ‘ x o _ y : yp _ z - z 0 _ t 

0 sin a cos a 
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is the equation of a solar ray parallel to the unit vector u, the interception 
of this straight line and the surface of the paraboloid is defined by the 
following s ?t of equations : 


x 


y 


z 


z 


X 


0 


y^ + t sin 

Zq + t cos a 

2 2 

M* + y )• 


(87) 


Here we have four equations with four parameters: Xq, y^, Zq, and t. 

To be able to eliminate them we must introduce the condition that the straight 

2 2 

line is tangent to the surface z = \(x + y ). Thus we have 


Zq + t cos a 



+ t sin 



or 


(X sin 2 a)t 2 


(2\y 0 sin a - cos <*)t + |x.|x 2 + Yq) " z o] = 0 ’ 


a quadratic equation in t which, due to the condition that the straight line 
has to be tangent to the paraboloid, has to have equal roots, t^ = t^. Since 


cos a 


1/2 


2Xy Q sin a ± ^cos a - 2XyQ sin - 4\ 


0 

2X sin^ a 


. 2 

sin a 


(*o * yp) - *o 


we have 


*1 = *2 = 4 = 


cos a - 2X.yQ sin a 
2 a . s m a 
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and 


(cos a - 2>.yQ sin a) = 4X. sin 


“ X ( x 0 + tS) - z o_ • 


( 88 ) 


The substitution of the above- given value for t into the second of Eqs. (87) 
yields 


cot a 6 . . , ^ 0 

y = = 2 cot cot or = — = y g , (89) 

i.e. , the value previously given by the second of Eqs. (44). This equation 
shows that the limit of the shadow is a plane curve, contained in a plane 
parallel to the z-axis, which projects in the xy- plane as the part of a straight 
line parallel to the x-axis. From the first and third of Eqs. (87) and from 
Eq. (88) we can derive the equation of the limiting curve in the xz- plane, 
eliminating Xq and Zq. Setting y = yg, we find that the curve is the parabola 

2 y S 

z = X.x + — cot or , 

i.e. , the parabola z = Xx^ (projection of the reflector on the xz-plane) moved 
translationally upward. The projection of the illuminated surface area of the 
reflector is shown on Fig. 12. 

The gradient of the outside surface of the reflector is pointed outwards; 
in other words, its direction is opposite to the direction of the gradient of the 
inside surface of the reflector. Thus, from Eq. (85), 

V& ~ 2\xe"j + 2\ye*2 - , 


and, from Eq. (86), 


N* = — (2\xe^ + 2Xye2 - e^) 


> 
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where W was given by Eq. (51). The cosine of the angle 0 between the normal 
to the surface $ and the direction of the incoming radiation is then 


z* - rT* 2\y sin a - cos ct 

cos 0 = u • N* = — 1 ^ . 


(90) 


From what was previously said, we orient the normal to the surface in 
such a manner that the angle 0 lies in the first quadrant; i. e. , so that 


cos 0 = u * N* > 0 


This condition yields 


2Xy sin ct - cos a > 0 . 

For 90° < a < 180°, y > 0, sin a > 0, - cos a > 0, and the condition is 
always satisfied. For 90° - ft < ct < 90°, the inequality above yields 
the constraint 


cot a 

y > ~i r = * 

which is always satisfied for the illuminated portion of the back surface of the 
reflector. 

The solar radiation force on the back surface can be found in the same 
way as was done for the front side of the reflector. For 90° - ft < a < 

90° + ft, one part of the outside surface is illuminated, while for 90° *t 
ft < a < 180°, the whole back side is lit by solar rays. The force over the 
whole area can be obtained from the corresponding expressions for the front 
side of the reflector if we note that cos 0 from Eq. (90) is the negative of the 
expression for cos 0 used for the front side, given by Eq. (53). Thus, denoting 
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by a prime the integrals corresponding to integrals I^, I^, 1 ^, I ^ and I 
for the front surface, we find 


I' = I 
1 21 21 


I* = -I 
l 21 l 22 


I' = I 
31 31 


(91) 


r = - 1 

a 32 32 


i = i 


and the components of the solar radiation force will be 


f'° = 0 

x 


F = K(-C,I-. + C.I,, - C-.I sin or ) 
y 321 1 22 2 


(92) 


F = K(C,I,, - C.L, - C_I cos a) 
z 3 31 1 32 2 


where 132 * and I are given by Eqs. (63), (64), (65), ( 66 ), and 

(67), respectively. 

Equations corresponding to Eqs. ( 68 ) for the front side are then 
F ^ |f^ sin o> + f 2 sin 2 aj 


iq 2 l ' ' ' \ 

F z = -it 6 K lg Q + gj cos a + g 7 cos 2ol 


(93) 
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where 


so that 



«<> = - g o 
g i = g l 
g 2 = " g 2 


F ^ = -ir6^K^ sin o - f^ sin2oj 

' 02 / \ ' 
F z = -it 6 K|-g Q + gj cos a - g 2 cos 201 


( 94 ) 


(95) 


where fj, f 2> g Q , gj, and g 2 are given by Eqs. (69), (70), (71), (72), and 
(73), respectively. 

Introducing again the ratio 


n 1 ( Q ') 


Projection of the illuminated area 
Area of the circle 


we find, for 90° < a < 90° + f2, from Fig. 12, 



2<(> 0 + sin 2<f> 0 \ 
¥ / 


(96) 
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so that the total force on the illuminated portion of the back side is 
approximately 

F y = 1lO%° 

■ n,«»F; 0 


( 97 ) 


where 


- 


'll 24» 0 +sin2<t> 0 \ 

|ll - — ^ -I for 90° - £2 < n < 90° 


+ £2 


1 for 90° < o < 180°. 


The graph of this function, for the Pioneer Spacecraft, is shown on Fig. 13. 
From Eqs. (78) and (96) we can also write 


i^) 



(98) 


This method, however, does not give as good results as in the case of the 
front side of the reflector and can possibly be used only for the values of a 
which are not too far below the value a = 90° + £2. For example, taking 
a = 90° and assuming black body structure of the back side of the reflector 
we find F' = F^ = 0, which is not true. Therefore, the components of the 
radiation force in the whole interval 90® - £2 < a < 90° + should be determined 
by the double integrations over the area of projection A j 


F' 


x 


= 0, 


F’ 

y 


F' 

z 




2 

sin $ dr d<j> - KC 2 I sin a , 


xy 


J5 

4\ 


jj B(6)r 


dr d<f> - KC 2 I cos a 


xy 
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■where B(6) is given by Eq. (21), cos 6 is given by Eq. (90), and 


I 




cos 0 r dr d() . 


The integration is simpler than in the case of the front side of the reflector 

because the projection of the illuminated area is the segment of the circle 
2 2 2 

x + y = 6 . The integration limits are 
For 4> = 4 >q to 4 > = tt - 4>q 

sin 4> 0 

For r: r = 2m — : — — to r = 2m , 
s 1119 

where <f> Q is given by Eq. (90). 
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Fig. 11. Solar radiation on the convex side of the 
parabolic reflector 
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^(a) 



Fig. 12. Projection of the illuminated outside area of the 
reflector on the xy-plane of reference 



Fig. 13. Graph of the blocking function for the convex side of the 
parabolic reflector 
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VIII. THE MOMENT OF THE SOLAR RADIATION FORCE ON THE 
BACK SIDE OF THE PARABOLIC REFLECTOR 


In the same manner as for the inside area of the reflector, we can 
prove that 


M' = M' = 0. 

y 2 


Since, for the back side, 


|N* • e 3 | 


= -(ye , -xe-)(l + 2\ z)dx dy f 


and 


uXr = Cj(z sin or - y cos a) + e^x cos a - e^x sin or, 


the expression for the moment, given by Eq. (31), becomes (note that for the 
back side r X dS is the negative of the corresponding value for the front side) 


M' 

x 


■//■ 


if 


= K II yB(0)(l + Z\z) dx dy + KC^ I I W(z sin a - y cos a) cos © dx dy, 


xy 


xy 


and, since cos 8 has the opposite sign of the same expression for the front 
side of the reflector, integrating over the whole area of the circle, we find 
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or 

M ® = ir6^CK(hj sin or - 1^ sin 2ar) (99) 


where hj and 1^ are given by Eqs. (83) 
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IX. TOTAL SOLAR RADIATION FORCE AND TORQUE ON BOTH SIDES OF 
THE PARABOLIC REFLECTOR: RECAPITULATION OF FORMULAS 


In the general case, the reflectivity characteristics of the front and 
back sides of the parabolic antenna reflector are different. We shall denote 
by Yp the values of p and y for the front side, and by pg, y fi the values 
of the same quantities for the back side of the reflector. Also, set 


py = p 


so that 


^F = P F Y F 
^B = P B Y B 


( 100 ) 


( 101 ) 


Let us compute first, for given values of 6 and ^ , 


ft = 


arctan 



and the constants of the reflector 


a^ = cot ft In cos ft 


a, = 


_ 2 1 - gpsJj (2 . ot 
“ 9 1 + cos ft + sec 


4 cos ft 


31+ cos ft 


b l " 15 


2 3 

2 1 + 2 cos ft + 3 cos ft + 4 cos ft 


cos ft(l + cos ft) 


b 2 = cot ft . 


( 102 ) 
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Then 


for the front side, and 



,B 

1 1 = a 2 B 

f 2 = ? + (? + 

B 1 

g 0 = 2 ' a l fl B 
B 

®1 = a 3 v B 

e? - < 1+3a i>>*B 

h i = Vb 
h 2 * \ + “iVb 



(103) 


(104) 
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for the back side of the reflector, where 


= Yi 


B 


Vd “ ^1 


( 105 ) 


The components of the solar radiation force and torque are schematically 
shown on Fig. H, ior different intervals of the angle a between 0 and 180 deg. 
If we redefine the functions T(o') and t)j(c) in the following manner 


1, for 0° < a < 90° - fl 


j 

={t f( 2 +n + sin for 90° -ft < a < 90 c 


0, for 90° < a < 180* 


we can write, for the components of the solar radiation force and torque, the 
following expressions: 


T 0 f \ 

f 1 = n(»)F + r f 

y y y 

T 0 1 

f x = n(tt)F + r f 

z z z 

M T = ri(o)M° + rM ' 
x x x 

where 

( 0 for 0° < a > 90° - n 

H 

1 for 90° - J2 < a < 180° 


(106) 
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Figure 14 shows the values of the components of the solar radiation 
force and moment for different values of angle a, as well as the values of 
the total force on both sides of the parabolic reflector. 

Equations (106), together with Eqs. (103), (104), and (107), represent 
the final form of the solution. 
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(o) FORCES AND TORQUES ON FRONT AND BACK SIDES 



FRONT SIDE 

BACK SIDE 

FRONT 

BACK 

Cl 

F 

F 

F 

F 

F 

F 

M 

M 


X 

y 

z 

X 

y 

z 

X 

X 

[o% 90° -a) 

0 

F° 

y 

F° 

z 

0 

0 

0 

M° 

X 

0 

(90° - Cl, 90°) 

0 

F 

y 

F z 

0 

F' 

y 

K 


< 

(90°. 90° + Cl) 

0 

0 

0 

0 

f ; 

F * 

0 

m' 

X 

[90° + fl, 180°] 

0 

0 

0 

0 

F '° 

F'° 

0 

m'° 






y 

z 


X 


(b) TOTAL FORCE AND TORQUE 


a 

F 

X 

F y 

F z 


0 
0 

$ 

1 

so 

0 

F° 

y 

F° 

z 

M? 

(90° - Cl, 90°) 

0 

F + F' 

y y 

F z + F z 

M * + < 

[90°, 90°+ ft) 

0 

f ; 

F z 

M / 

X 

[90® + n, 180®] 

0 

f'° 

F /0 

M /0 



y 

z 

X 


Fig. 14. Table of the total radiation force and torques 
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X. COMPARISON WITH THE FLAT- SURFACE MODEL 


The so-called "effective area" of the illumination in the flat surface 
model is the projection of the illuminated area on the plane perpendicular to 
the direction of the solar radiation. To simplify the derivations of the flat 
model (the model used in Ref. 5) and stiU be able to find the order of magni- 
tude of the deviation of the flat surface model from the real force model, we 
shall consider the case when 0 < a S 90° - f2. We shall also assume that 
the values of (3 and y are the same for both models. Then the components of 
the solar radiation force for a flat surface can be derived from Eqs. (68) by 
setting £ = 0. Then £2 = 0 and 


fj = 0, 


U = + 2py lim | > ° T- a - 

a— o[ 


-x i L~yj y liiJU I ^ “ 

~ " 1 tan , 


1 - p v 

~2 — ’ 


gn = 


l±M 

2 


so that 


^ - —%• (1 “ 3), 


1+gV 

g 2 2: 


0 2 

(flat) = -tt 6 K(1 - p-y) sin a cos a , 


F^ (flat) = -it6^k£^( 1 - p) cos a + (1 + Pv) cos^ arj- 


The differences in the components of the real force and the flat surface 
force model are 
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= -*5 2 k[><. (1 . p) } ; fgfg (2 + seen) sin or 

+ P"y ( 1 + 2 cot 2 nin cos n) sin 2a J 

A ^0 c 2 t ^P2y /i 0 , 1 - cos n 
AF z = *6 (1 - (3) T - CQS n cos « 

+ ^ (1 + 3 cos 2a)(l + 2 cot 2 J2 In cos n)] 


For the Pioneer F/G spacecraft, these differences are 
AF° = -tt6 2 K [o. 046754 y(l - p) sin a + 0. 128078 Py sin 2aj 

AF® = *6 2 k[o. 044621 y(l - P) cos a + 0.064039 Py( 1 + 3 cos 2a)J 

which shows the order of magnitude of the error in the flat surface model. 
For example, for <* = 0, p=*y = l, 

AF° = 0. 256156 tt6 2 K 
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XI. PRINCIPAL AXES OF INERTIA; EULER'S EQUATIONS OF MOTION 


In order to derive the equations of motion of the spacecraft under the 

action of the solar radiation force, we shall assume that the spacecraft is 

rotationally symmetric and that its center of gravity (point C on Fig. 15) lies 

on the z-axis. Now define the system of principal axes of inertia Xp, yp, Zp 

of the body in the following manner. Because of symmetry the principal 

z -axis will coincide with the z-axis of one frame of reference and, for the 
P 

same reason, the x^ and y^ axes may be chosen arbitrarily. Therefore, to 
simplify he problem, we shall take the x^ and y^ axes to be parallel to the 
x and y axes respectively. Hence, 

x = x p 
y = yp 

Z = - Z^. 


The components of the moment of the solar radiation force with respect to 
the point O, the origin of the system xyz, are 


M x ( °> = M x 


My 


(O) _ 


= o 


M z {0) = O. 


If r is the position vector of a point of the body with respect to O as the 
—(C) 

origin and r is the position vector of the same point with respect to the 
point C as the origin, the relationship between these two vectors is 


r^) = r + z e-> . 

C 5 

The moments of the solar radiation force relative to the points O and C are, 
respectively. 
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(o) 


II r X 

s 



fjr^ X dF = 

S 



+ Z c e 3 


xF. 


further on, because F = O, 

x 


e, X F = -F e, 

3 y 1 

so that 

M (C * = M - z „F e, 
O y 1 


or, consequently, 


M < G > . M 

X X 


M 

Y 


(C) = 


O 


M 

z 


(C) 


o 


z r F 

c y 


where 


x = x 
P 


y 


p 


y 


z 

p 


z 


z 


C 


are the principal axes. 


(108) 
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The moments of inertia of the hollow paraboloid of revolution with 
respect to the axes of the system xyz are 


I 


x 


I 

y 



+ 3 ; 2 ) 


I 


Z 



where is the mass of the reflector. According to Steiner's theorem, the 
principal moments of inertia of the reflector are 


A = 



(109) 


Let J , J , J be the principal moments of inertia of the rest of the 
xyz r 

spacecraft. Then, using Eqs. (108) and (109), we can v/rite Euler 's equations 
of motion of the whole spacecraft in the form 


(A + J x ) ^ - [(A - C) f (J y - J z )] * y * z 
(A + J y ) iy + [(A - C) + (J x - J s )] % 



(HO) 


Here we assume that the rest of the spacecraft is not illuminated (w , ui , w 

x y z 

are components of the angular velocity vector). 

For a rotationally symmetric spacecraft, = J , and 


to 


z 


n 
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so that 


A + J 

CO = - Gj = - £ (Jj 

x n[( A - o + (j x - j z >] y y 


and the differential equation for u)^ is 


** 1 J ** ii + n (A - C + J - J ) u = z F - M. (Ill) 

„(A-C + J - J ) y x z' y C y * 

X Z 


The right-hand side of this equation can be written in another form, by means 
of the first of Eqs. (68) and Eq. (82), namely 


z- F - M 
C y x 


= - tt 6 2 K [(fj z c + 4hj) sin a + (f 2 z Q + Ch 2 ) sin 2aJ. (112) 


From Eqs. (69), (70) and the first of Eqs. (83), we can find the values 
of f j z q + £hj and f^ z ^ + £h^ in terms of the angle S2 of the reflector. 

F rom 


, m k srp 
K - K ( r gp) = 


using the notation 


f i z c +Ch i = e i 


f 2 2 C + ^ h 2 ” 1 2 


(113) 


K SRP " C S 


we have the right-hand side of Eq. (Ill) in the form 


F — - -Cg (( j sin at t sin 2 or). 


(114) 
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For the Pioneer spacecraft, 

| s Y (1 . p) (0. 127367 + 0.046754 z Q ) meters 

l z = i (0. 3803 + z c ) + (0. 158396 - 0. 371922 z Q ) PY meters 

17 2 

C_ = 6.0444 X 10 1 ' kilogram meter /second 

O 

or, if the heliocentric distance of the spacecraft is given in astronomical 
units (AU), 

Cg = 2.7008 X 10“ ^ kilogram meter /second^ 

With a negligible loss of accuracy we can assume that the angle a is equal to 
the Earth-Spacecraft-Sun angle. Therefore, 

|* XR E | 
sin « = — 

Using the notation 

n (A - C + J x - J z ) = q j 
(A ♦ j/ = q 2 

we can rewrite Eq. (Ill) in the form 
92 C S 

w + q, w = - (£, sin a + £ sin Za) 

y i y r 2 1 2 


or 


W + 

y 0 


u = 


€^(1 + e cos f) 


(I, 


sin a + #2 sin Za) 


(115) 
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where 


% ' 


n (A - C + J x - J z ) 
A + J 


0 ag (1 - e 2 ) 2 (A + J x ) 2 

= semimajor axis of the spacecraft's orbit in astronomical units 
e = the eccentricity of the spacecraft's orbit 
f = true anomaly of the spacecraft- 

The period of free oscillations (without the forcing function) of the rotationally 

symmetric spacecraft (J = J ) is given by 

x v 


c s n 


A - C + J x - J z 


, A + J x 

T _ 2tt x 

1 0 ” n A-C+J-J 

X z 


(116) 


where n is the rate of rotation of the spacecraft about its roll axis. 

If there were no other parts of the spacecraft (the reflector only). 


J = J = 0 
x z 


Z C = 3 * 


A M (£,± ±£\ M (§L + sZ\ 

A-M R\3 2"9y y 3 + 18/ 


M R T 


A - C = 


i Z 

M R k 
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and 



For the Pioneer spacecraft, the rate of rotation about its roll axis is 
5 revolutions/minute, and the period of spinning is thus 

T = 12 seconds, 
spin 


Therefore, 


T 0 = T spin [' ♦ 6 (f) 2 ]- 


The ratio 6/£ is approximately 3.6, so that 


T_ = 79 T . = 948 seconds = 15.8 minutes. 

0 spin 

The period of free oscillations about the x-axis is the same. Indeed, setting 
the right-hand side of the first of Eqs. (110) equal to zero, we can substitute 
into the second of Eqs. (110) and obtain again Eq. (115) in the form 

w + flj u =0. 

x Ox 
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XU. SOLAR RADIATION FORCE IN THE INERTIAL FRAME OF REFERENCE 

The components of the solar radiation force, given by the first two of 
Eqs. (106), are relative to the noninertial frame of reference xyz. The 
angular velocity vector of rotation of this system is given by 

2w = Tx-gJ +T x -^+ kX^ (117) 


where 


X r 


i = 


I b e x7 I 


j = kx i 


k = 


Here, 




‘ x dt *T- S ~7*'h (r e x r> 

dt I*e x7 | dl * 


— dT X (®e* xr) j p . 

' " (p|r_x 7|) 2 X * K ' r ’ X ( * eX r ’] 


j X 


r- x & R E- r x /dR 

kX X = - 72 _X 


' dt dt / 


and the equation of motion of the spacecraft in an inertial frame of reference 
is 


.2— _ j— _ _ J_ _ Tp 

■ a ~r- - 2u X ~ + u X (u X r) - w X r = — — 
dt 2 dt m S 


where F = F^j + F^lc, and mg is the mass of the spacecraft. 
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XIII. SOLAR RADIATION FORCE ON THE SURFACE 
OF A CIRCULAR CYLINDER 


A spacecraft may carry one or more components of cylindrical shape 
such as fuel tanks, instruments, etc. To find the solar radiation force 
exerted on the surface of a cylinder of an arbitrary radius in the general 
case, we shall assume that the axis of symmetry of the cylinder lies along 
an arbitrarily chosen direction, defined by the unit vector IT ^ . Let xyz be 
the noninertial frame of reference, the axes of which are directed along the 
unit vectors e“j , "e^, ^"3 the already described system of coordirate^, trans- 
lated into the center of gravity as the origin, and let the point O be the ^oint 
on which the axis of symmetry of the cylinder penetrates the yz-plane of the 
above -described coordinate system. Taking O as the origin of another coor- 
dinate system of axes x', y*, z‘ parallel to the axes x, y, and z respectively, 
we have, from Fig. 16, 


x' = x 

y' = y - y 0 


where the two quantities and z^, are known# In order to bring the axes of 
the system x* y 1 z' into the position of the system £ r) £ , we must perform the 
following two rotations: 

(1) Rotation about the z'-axis by an angle y, in the positive direction; 
this brings the x'-axis into the position OL. 

(2) Rotation about the new y-axis by an angle € , in the negative direc- 
tion. This brings the x*-axis into the position 0£. 
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The equations of transformation between the x' y' z 1 system and £ q £, 
system are 


" 5 ' 


COS € 

0 

sin c 


cos x 

sin x 

0 


x'" 

q 

= 

0 

1 

0 


-sin x 

cos x 

0 


y’ 

4 


-sin c 

0 

COS € 


0 

0 

1 


z' 


or 


4 


COS X COS € 

sin x cos € 

sin € 


x' 

q 

= 

-sin X 

cos X 

0 


y' 

C 


-cos x sin € 

-sin x sin « 

COS € 


z' 


(119) 


The equation of the cylinder in the system | q £ is 


2 , ,2 
n + ; 


(120) 


The components of the unit vector 

u = e 2 sin a + cos a 
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in the |r)£ -system are 


u, 


U, 


U. 


cos x cos € sm x cos < sin c 


-sin x 


cos x 


-cos x sin « -sin x sin c cos c 


sin a 


cos a 


or 


u, = sin x cos < sin a + sin < cos or 


= cos x sin O' 


u, = -sin x sin c sin a + cos c cos a. 


and the equation of the plane parallel to unit vector u and n^ is 


(r - r Q ) • (uXBj) = 0 


or 


£ " So 1 ” To 6 ' 4 0 


U 1 u 2 u 3 


= 0 , 


( 121 ) 
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or 


u 3 (n - n 0 ) ■ u 2 ^ = °- 

This equation can be written also in the form 

u 3 q - ^ + p = 0. (122) 

Out of infinitely many planes given by Eq. (122), two are tangent to 
the surface of the cylinder. Their traces are shown on Fig. 17. The two 
lines along which these planes touch the cylinder project a two points. P 
and Q, on Fig. 17. These two lines separate the illuminated part of the sur- 
face from the part of the surface of the cylinder which remains in the shadow. 
The illuminated area projects on the ^rj -plane as a rectangle of length AB = h 
and width 2rj ^ (Fig. 17), whiv-h must be determined in the following way. 

Combining the equations of the cylinder (Eq. 120) with the equation of 
the plane (Eq. 122), we find the interception of these two surfaces to be given 
by 


n 



. JL 
u 3 


Substituting the value of q from the first equation into the second, we obtain the 
the quadratic equation for C: 



0 . 


Because of the condition for tangency, the two roots of this equation must be 
the same. Hence, we find 


P 


* a + u 3 



(123) 
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because 


2 , 2 , 2 . 
u 1 + u 2 + u 3 = 1. 


Therefore, the coordinates of points P and Q are 


au., 

S>. Q ^ 


V 1 ' u i 


and 


n - ± 
P, Q 


au. 




so that 


T = 
0 


au. 


yf 7 * 


(124) 


The unit vector of the normal to the surface 


— V® 

N *=iv*i =i r 


lies in the q £ - plane. From the equation of the surface 


«(t. n, C) = *i 2 + C 2 - a 2 = 0 


we find 


3 ® 9 ® 9 * 

if = °- an " 2l1 -^ * 
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so that 




W = \/40l 2 + ; 2 ) = 2a 


and 


— n — — 

N* = -n, + — n 0 . 
a ^ a 3 


Also, from cos 6 = u • N*, we find 


cos 0 = 


u 2 ri + u 3^ 


Taking the coordinate transformation 

H = a sin ^ 
£ = a cos <t> 

we -an write, instead of Eq. (127), 


cos ® = u 2 S * T) + u 2 cos 


Since the H -coordinate of the point Q (Fig. 17) is 

n Q = n o * a sin A 

the integration limits for the solar radiation force in Eq. (29) will be 


for | : 

^A’ k 

for <)> : 

: -A, +A, 


(125) 


(126) 


(127) 


(128) 
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where 



The solar radiation force is then given by 


( 129 ) 


,.„r f*- h 

J(br,-A Jt=i 


■ a, _l — A\ a cos ^ d| d<j> 

B( 0 )(n ? sin 9 + n- cos o) — — 


_ _ _ l " '' r «A’ h a cos 4> d£ d 4* 

n + u 2 n + u 3 n ) / / cos 6 . jf,| ’ 


-KC 2 (u 


or, because 


|n 3 • N*j = - = cos 


■£ 


F‘ = -Kha/ B( 0 )(n 2 sin * + u 3 cos <|>) d<|> 


-KhC 2 a(u 


l n j + u 2 n 2 + u 3 n 3 )jf 


cos 0 d«j> 


( 130 ) 


where 


5 ( 0 ) - Cj(it 2 sin <|) + u 3 cog <fr) + 2 C 3 (u 2 


sin 4 > + u 3 cos 


4 >) 
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The components of the solar radiation force, given by Eq. (130), 
the axes of the system are 


'■■I 


= -KhaC^Uj I (u^ sin 4> + cos <J>) d6 

'-A 


( 


u. 


= -Kha j B(6) sin 4> d4> + ~ F^ 


u. 


F, = -Kha / B(6) cos 4> d<t> + — F . 

U 1 5 


Performing the indicated integration, we finally obtain 


F^ = -ZKahC^UjU^ sin A = -ZKahC^ 


U 1 U 3 




F^ = -Khau 


8C, 





F^ = -Khau ^ 


C ^ arc sin 




8C. 


i 2 2v 
1 - u l +u 2 \ 

V? I 


U 3 

+ — F, 

U 1 « 


along 


(131) 


(132) 


(133) 
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The components of the solar radiation force along the axes of the system 
x'y'z ' will then be 


~ 


“■ 




— 


— • 1 

F 

X 


cos X 

COS € 

- sin X 

- cos X 

sin € 


F t 

F 

y 

= 

sin X 

COS € 

cos X 

- sin X 

sin € 


F n 

F 

z 


sin 

€ 

0 

cos € 










__ 

L J 


(134) 


Since the system x'y'z' is rotating in space, and the spacecraft also spins 
about its roll axis, the angle X is a function of time. 
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Fig. 16. Solar radiation on the curved surface of a 
circular cylinder 



Fig. 17. Trace of the projection of the illuminated 
area of the cylinder on the xy-plane 
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XIV. SPECIAL, CASES OF THE SOLAR RADIATION 
FORCE ON A CIRCULAR CYLINDER 


If the circular cylinder lies in the x'y' plane, 6 = 0, and, from 
Eqs. (121), 


Uj = sin X sin a 


U£ = cos x sin or 


- cos a 


so that, from Eq. (129), 


sin A 


cos or 



. 2 

sm X 


. 2 

sin a 


From Eq. (131) we can derive 


3 

sin x sin a cos a 

/ 2 2 
V 1 - sin X sin a 

and the other two components can be obtained by a simple substitution of the 
values of Uj, u^, and u^ into Eqs. (132) and (133). 

Another interesting case is noted when the axis of symmetry of the 
cylinder coincides with the z-axis (axis of symmetry of the reflector). Then, 
the |-axis lies along the z-axis, the q-axis coincides with the x-axis, and 
the 1,-axis coincides with the y-axis, so that X = " 7T /2> and € = tt/ 2. Assum- 
ing that the cylinder hangs behind the surface of the parabolic reflector, we 
can take the bottom base of the cylinder as the q£ -plane, which yields 
y^ = 0, Zq = h. Thus 
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' U 1 


0 0 1 


0 

u 2 

= 

1 0 0 


sin a 

. U 3 


0 1 0 


cos a 


or 


Uj = cos a 



= sin ot. 


Neglecting the shadowing effects of other components of the spacecraft on the 
surface of the cylinder, we can write from Eq. (129) 

sin A = 1, A = 90°. 

Formulas given by Eqs. (131-133) yield 

F^ = -Kah(l - py) sin 2a 

F n = 0 

= -Kah sin ap^(l - P) + •j(3 + p-y) sin orj 
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which leads to the following expressions 



F = 0 \ 

x 1 

Fy = -Kah sin (1 - P) + y(3 + Py) sin arj \ (135) 

F = -Kah(l - Py) sin la. ) 

z 
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APPENDIX 

STOKES' THEOREM; CONVERSION OF DOUBLE 
INTEGRALS INTO LINE INTEGRALS 


The previously obtained expression for the solar radiation force, given 
by Eq. (22), may be written in the form 


F = 


-K 


// 


B(e)dS + C 2 u 



dS 


(A- 1) 


i. e- , as the sum of two double integrals. We have seen earlier that the 
second integral in Eq. (A-l) can be solved immediately, since, insofar as 
the double integration is concerned, the unit vector u is a constant vector 
and, therefore. 



(A-2) 


The function B(0) in the first double integral in Eq. (A-l) is a function 
of x, y, and z. From the equation of the reflecting surface, 

$(x,y,z) = 0 


we can obtain z as a function of x and y. With the substitution of the value of 
z so obtained into the function B(0), it becomes a function of x and y only, 
i.e. , B(x,y). Thus we have 


F 


-K 


'// 


B(x,y) dS + C£(u • S)u 


(A- 3) 
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Now consider the double integral 


B(x, y) d3 . (A-4) 

Since this integral physically represents one part of the solar radiation force 
exerted on the surface of a certain component of a spacecraft, the surface 
area S is always bounded by a simple closed curve or, in the general case, a 
closed geometrical figure bounded by several parts of different simple 
curves. Therefore, we can write the Stokes theorem in the form 



if 


(dS X v) X v = 




X dr. 


(A- 5) 


where L is the closed curve, bounding the surface S. The vector v is an 
unknown vector function that must be determined from 

(dS X V) X v = B(x, y) dS . (A-6) 

Since, from what we had before. 


— VCD 

dS = N* dS = dS, 


we can rewrite Eq. (A-6) in the form 

(V$ X V) X v = B(x,y)V$. (A-7) 

The equation of the surface $ having been known, we can consider a known 
vector , 

N = 
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the components of which are N^(x,y), i = 1,2,3. The left-hand side of 
Eq. (A-7) is a triple cross-product, which can be written as 

* x * 

(N X V) X v = V(v • N) - N(V • v) = grad(v • N) - N div v (A- 8) 


where the asterisk abcve the vector N means that the operator V is applied 
to vector v only. 

Substituting the value of the triple cross-product (N X V) X v from 
Eq. (A- 8) into Eq. (A-7), we find 

* 

grad(v * N) - N div v = B(x,y)N 


or 


grad (v • N) n N^iiv v + B(x,y)J 


(A-9) 


Using components of the vector function v: v^, v^, and v^ t Eq. (A-9) yields 

the following three partial differential equations: 


9v. 


dv , 


+N 2^T +N 


fn . N fci 

3 9x 1 L 


9v ? 9v^ 

+ TT + aT + E|x 


3v, 9v ? 9v^ 

N i af + N 2a7 +N 3a7 


9v, fav, 3 v 2 9v, *| 

W - N 2L3F + W + +B <*’ y 'J 


(A- 10) 


9v, 9v ? 9v- 9v. dv? 9v~ 

N i aT +n 2-5T +n 3 3F * N jb* + W + + B<x ’ y) . 


or 
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[*3v ? 3v_ *] 3v ? 3v, 

N ila7 + ?r + *<*•*>} - n 2 3T - N 3ir * 0 


9v i f 

- N 1 37 +N 2L 


9v l ^ v 3 
3x + 3z + 


1 9v 3 

b(x» y)J - N 3 'a7 = 0 


3v. 3 v 2 T® v 1 ^ v 2 

N 1 -5T - n 2 3T +N 3Lar + W + B(X 


,y)l = o 


where B(x,y) is a known function of x and y. Since this system is homogeneous 
with respect to N^, and Nj , the determinant of the system must be zero. 

We shall now consider a particular case when the right-hand side of 
Eq. (A-9) is zero; in other words, we shall find a particular vector function 
v which makes the right-hand side of the Eq. (A-9) vanish, so that 


div v + B(x, y) = 0. 


lA-11) 


Since in this case the right-hand sides of the three scalar equations 

(Eqs. A- 10) vanish, the Jacobian J(v^, v^, v^/x,y,z) is equal to zero, i.e., 


J 


/ v r v z» v 3 \ 

V x,y,z / 


d Zi 

d Z 2 

^3 

9x 

9x 

9x 


^2 

^3 

3y 

3y 

3y 



3V3 

3z 

9z 

9z 


This fact implies that the three components of the vector function v are not 
independent quantities. Indeed, writing the equation of the reflecting surface 
in the explicit form 


*(x,y, z) = f(x, y) - z = 0 
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and, therefore, 



We find, from the third of Eqs. (A- 10), 


_9 /9f Of 

9z\3x V l 9y V 2 



0 . 


and, accordingly, 

3f X af , 

9x V 1 + a 7 v 2 - v 3 = 

or 

v 3 = ¥x v i + § 7 v 2 • (A ' 12) 

Hence 


3v 3 

a 2 f 

+ 2jL v + 

gi!ii 

9f^ 

d'P 

3x ‘ 

= ax 2 " 1 

dxdy 2 

3x 9x 

ay ax 

dx 

9v 3 

. j£l 

x ^ + 

t»h , 

. Sitli . 

d± 

dy ‘ 

" 3xay v 

1 » 2 2 
»y 

3x 3y 

dy dy 

dy 

»3 . 


* 3f ® V 2 




3z 

= 

3x 3z 

3y 9z 
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and, substituting these values into the left-hand sides of the first two of 
Eqs, (A , we find 


9 Z f , a 2 f 

9x Z v l 8x9y v 2 


9 2 f + 9 2 f 

3^ Vl + ^ 



(A- 13) 


Assuming that the determinant of this system does not vanish, i.e. , if 


G(x.y) = 


it 

9x 2 


9 Z f 

9x9y 


A 

9x9y 


9 2 f 


9y 




dx dy 


the two components v^ and v 0 can be determined from 


v, = 


1 / 9 2 f 94- 

G(x,y) \9^ dx 


9 2 f 94- 
9x9y 9y 


) 


1 / 9 2 f 94- 9 2 £ 94 \ 

G(x,y) 9y ' 9x9y9x) 


(A- 14) 


Eqs. (A- 14) show that Vj and v^ are functions of x and y only. Therefore 


9vj ^ 9v^ 
9z - 9z 


= 0 
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and, accordingly, 



The conditional scalar equation (Eq. A- 11) then yields 



dx 



-B(x, y). 


(A- 15) 


so that, differentiating the first equation (Eq. A- 14) with respect to x and the 
second equation with respect to y, and substituting the results into Eq. (A- 15), 
we finally arrive at a partial differential equation of the second order lor the 
determination of the function (x,y). 

To illustrate this method, we shall consider the equation of the surface 
of the paraboloid of revolution 

$(x, y, z) = Mx 2 + y 2 ) - z = 0 


where 


8f 3£ 

- = 2\x,^ = 2Xy 


a 2 f ,, a 2 f 

ax 2 = ’ 9xa y 


o. 


a 2 f 


= z\ 


G(x, y) = 4X. 


Here, from Eqs. (A-14), 


1 aji, 

V 1 " 2 \ 9x 


) 

v 2 " 2T 8y 
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and Eq. (A- 15) yields the Laplace equation 


Aip 


iLi + jLi. _ 

dx 2 dy 2 


in two dimensions: 
-2X.B(x, y). 
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NOMENCLATURE 


a 


a 0 


a r a 2* 


a 3 


a 


1 * 



semidiameter of the circular cylinder 

semin. i.jor axis of the spacecraft's orbit 

constants oi the parabolic reflector 

unit vectors along the axes of the inertial reference 
frame 


a SRP 

A 


acceleration of the solar radiation force 


moment of inertia of the parabolic reflector about 
the principal x^-axis of the spacecraft 

parameter used in integration 




projections of the illuminated surface area on the 
three coordinate planes 


A(f) 


reflectivity function 


AU 


astronomical unit 


B 

B(e) 


B(f) 

c 


C 



constants of the parabolic reflector 

moment of inertia of the reflector about the princi- 
pal Yp-axis of the spacecraft 

function of the angle between the direction of the 
incident radiation and normal to the reflecting 
surface 


reflectivity function 
speed of light 
modified solar constant 


moment of inertia of the parabolic reflector about 
the principal z^-axis of the spacecraft 


reflectivity characteristics of a specified reflecting 
surface 


e 


eccentricity of the spacecraft's orbit 


e 


r 


V 


e 3 


unit vectors along the axes of the noninertial frame 
of reference xyz 


f 


true anomaly cf the spacecraft in its orbit 
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2 

rF 



f(e) 

F 


constants of the parabolic reflector 

values of constants f^ and f^ for the front side ot 
the reflector 

values of constants f^ and f^ for the back side of 
the reflector 

reflection law 

solar radiation force 



F 


N 


F 


I 


F 


R 


F 


D 






components of the solar radiation force along the 
axes of the noninertial frame of reference xyz 

solar radiation force exerted upon a surface per- 
pendicular to the direction of the incoming 
radiation 

part of the radiation force due to the incoming 
radiation 

part of the radiation force due to the specularly 
reflected radiation 

part of the radiation force due tu the diffusely 
reflected radiation 

solar radiation force exerted on the whole front 
surface of the parabolic reflector 

components of force F^ along the axes of the non- 
inertial reference frame xyz 

components of the solar radiation force exerted 
on the whole back area of the reflector, along the 
axes of the noninertial reference frame xyz 

components of the total solar radiation force in 
the system xyz 



F n’ 


components of the solar radiation force along the 
axes of the noninertial system § d £ 




constants of the parabolic reflector 


F F F 

go* 8 1 g 2 


values of constants g~, g,, and g~ for the front side 
of the reflector 


BBT. 
’ g l * g 2 


values of constant gg, g^, and g ^ for the back side 
of the reflector ^ 
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h 



i 

I 

*0 

If *21* l zv hi* hz 


i 

z 


j 

J 


J ’ J 21’ J 22’ J 31’ J 32 


J , J , J 
x y z 


K = K(r gp ) = K(r) 


Kj, K 2 , K 3 


K 


SRP 


V l z 


m 


m r 


length of the circular cylinder 

constants of the reflector; coefficients in the 
expression for the moment of the radiation force 

values of constants h^ and h^ for the front side of 
the reflector 

values of h^ and h^ for the back side of the 
reflector 

unit vector along the x-axis 

radiant flux per unit solid angle on a hemisphere 

constant, having the same dimensions as I 

definite integrals that appear in the expressions 
for the components of the solar radiation force 

moments of inertia of the reflector about the axes 
of the coordinate frame xyz 

unit vector along the y-axis 

radiant energy of the incoming radiation per unit 
area per unit of time 

definite integrals that appear in the expressions 
for the components of the solar radiation force 

solar constant 

principal moments of inertia of the spacecraft 
without the parabolic reflector 

unit vector along the z-axis, in the Spacecraft- 
Earth direction 

radiation constant 

function of the heliocentric distance of the 
spacecraft 

reflecting characteristics of a surface 

solar radiation constant 

auxiliary constants 

constant of the reflector 

total mass of the spacecraft 
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M 


mass 


M 


moment of the solar radiation force 


M , M , M 
xyz 


M (°) M (°) 

M<°> 

x y 

z 

M (C \ M (C) , 

M< C > 

x y 

M°, M°, M° 
xyz 

z 


components of the vector M along the axes of the 
noninertial system xyz 

same as above; O is the origin of the system xyz 


components of the moment of the solar radiation 
force along axes parallel to x, y, z-axes, centered 
at the center of gravity of the spacecraft 

values of components M x> M , M z respectively 
for the whole inside area of tne reflector 


M r , M r , M' 
xyz 


components of the moment of the solar radiation 
force for the back side of the reflector 



M 

y 



values of components My, M' z for the whole 

back side area of the reflector 


m r 


n 


n l’ n 2’ n 3 

N* 


mass of the parabolic reflector 

rate of rotation of the spacecraft about its roll 
axis (z-axis) 

unit vectors along the axes of the system £ r\ £ 

unit vector along the normal to the reflecting 
surface 


P 


I 


component of force Fj along the normal to the 
reflecting surface 


P 


R 


component of force along the normal to the 
surface 


q r q 2 


r 




S 


S 


constants 

heliocentric position vector of the spacecraft, also 
r SP |r| r SP 
dimensionless polar coordinate 

position vector of a point of the spacecraft relative 
to its center of gravity 

heliocentric position vector of the Earth 

oriented surface 

surface area 
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T 


0 


period of free oscillations of the spacecraft about 
the x- and y-axis 


u 


unit vector in the Spacecraft- Sun direction 


u l’ u 2’ u 3 


components of the unit vector u along the axes of 
the coordinate system |t]£, 


W 


magnitude of the gradient vector 


x, y, z 
x', y', z' 


z 

P 


X, Y, Z 



a 0 

3 

Y 

6 


S, 

n 

n(a), ^(or) 

e 


noninertial, spacecraft- fixed coordinates 


principal axes of inertia of the spacecraft 

heliocentric inertial coordinates of the spacecraft 

heliocentric inertial coordinates of the Earth 

angle between the axis of symmetry of the reflector 
and the direction of incoming radiation 

Earth- Space c raft- Sun angle 

portion of photons reflected specularly 

portion of reflected photons 

semidiameter of the parabolic reflector 

angle of inclination of the axis of symmetry of the 
circular cylinder to the xy-plane 

depth of the parabolic reflector. Also, coordinate 
in the £r|£ system of reference 

coordinate in the£r|£ system of reference 

blockage (shadowing) functions for the front and 
back surface of the reflector, respectively 

angle between the normal to the surface and the 
direction of the incident radiation 




constat of the parabolic reflector 


A 


polar angle used in integration 




reflectivity characteristics of a surface 




values of p for the front and back surface of the 
reilecior, respectively 
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reflectivity characteristics of a surface 

values of v for the front and back surface of the 
reflector, respectively 

coordinate in the reference frame 

geocentric position vector of the spacecraft 



polar coordinate 

equation of the reflecting surface 

angle between the projection of the axis of sym- 
metry of the circular cylinder on the xy-plane and 
the x- axis 

angular dimension of the parabolic reflector 
angular velocity vector 
angular velocity 

components o'" the angular velocity vector along 
the axes of the noninertial frame of reference xyz 

angular dimension of the parabolic reflector 

frequency of the free oscillations of the spacecraft 
about the x- and y-axis 
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